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Abstract: 

Transport simulations including in-medium properties derived in a microscopic 7T+ NN^ 1 + 
A7V -1 model in infinite nuclear matter are presented. In-medium pion dispersion relations, 
partial A decay widths, pion absorption cross sections and A cross sections are incorporated 
into the transport description by means of a local-density approximation. Strong modifica- 
tions of 7r and A production and absorption rates are found, but only small effects on pion 
observables. 
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1 Introduction 



In collisions between two heavy nuclei at bombarding energies from a few hundred 
MeV up to several GeV per nucleon, hadronic matter at high density and temperature 
is formed and a large number of energetic particles are produced [|1], 0, |3|, ||| . The main 
features in such collisions have been fairly well described by microscopic transport 
models employing vacuum properties of resonances and mesons, such as most versions 
of BUU and QMD J2L [3], || . However, ir mesons, nucleons, and A isobars are strongly 
interacting particles, and in infinite nuclear matter a system of such particles would 
couple to form spin-isospin modes. It is quite possible that such in-medium effects 
could play an important role for the transport properties. First attempts to employ 
such in-medium modifications in transport simulations of nuclear collisions have been 
carried out. In Refs. JF], [8], |S], 110, |TT| the employed in- medium properties were based 
on a simple two-level AN' 1 model. A more complete tt + NN' 1 + AiV -1 model 
was used in Refs. |T2] to derive a consistent set of in-medium quantities suitable for 
transport descriptions, with exploratory transport simulations reported in Refs. |L3|. 

The previous attempts to include tt and A in-medium properties have suggested 
that only small effects are found in pion observables. This is not surprising since 
most emitted pions are created at the surface where the in-medium effects are small. 
However, in-medium modifications of pions and A isobars may be important for other 
observables. For example, presently large efforts are made to study effects of possible 
chiral restoration at high densities and temperatures, and in particular by investigat- 
ing kaon and dilepton spectra which may carry signatures of in-medium modifications 
of vector meson masses due to chiral restoration |n|, W , IB \. Since kaons and dilep- 



tons are produced through multistep processes, involving pions and A's, in-medium 
modifications of tt and A quantities could implicitly have an important impact. 

In this letter we present selected results from transport simulations that include 
in-medium quantities obtained from the tt + NN' 1 + AiV -1 model of Ref. The 



present work constitutes an extension of Ref. [13], now including a more complete 
and consistent set of in-medium quantities. Our intention is to present a qualitative 
picture of some in-medium effects that survive the transport dynamics, while a more 
complete and systematic presentation will be given in a subsequent paper. While 
the model presented here contains the main in-medium modifications of pions and 
A isobars, a number of additional processes might play a role as well. For example, 
effects originating from a possible partial chiral restoration have not been included. 
Moreover, there might also be in-medium modifications of cross sections and decay 
widths originating from the fact that the time and volume available in collisions and 
decays are finite. We expect that the present transport simulations contain the most 
important in-medium features. Furthermore, they provide the most consistent tests 
carried out so far of in-medium effects within AiV -1 models. 
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2 The Model 



Our approach is to obtain the medium modifications by microscopic calculations 
in uniform matter at various densities and then incorporate those into a transport 
treatment by a local-density approximation f[2 , For simplicity, all microscopic 
quantities are approximated by their zero-temperature values and the A width is 
omitted in the microscopic treatment of the dispersion relations. Detailed discussions 
of these approximations are given in Refs. |T2|, [T3[ . 

We consider interacting nucleons (N), delta isobars (A), pi mesons (it), and rho 
mesons (p) in a periodic box for T = 0. The in-medium properties are obtained by 
using the Greens function technique, starting from non-interacting hadrons. At pion 
vertices we use effective p-wave interactions containing form factors, F NnN (q) and 
Fnka(q)- Since various choices are made in the literature we will in this work utilize 
two different sets. The first set, denoted FF1, is taken 
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while for the second set, FF2, the relativistic corrections are replaced by an off-shell 
correction at vertices including a A isobar and neglected at other vertices. In addition 
the monopole form is approximated by an exponential form. Thus for FF2 we take 
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For most quantities of interest here the difference between the monopole and exponen- 
tial form is practically negligible, but the exponential form yields better convergence 
properties for calculating the real part of the A self energy. The relativistic cor- 
rection 2m a /[y/s + m a ] is frequently used in the A-hole model in connection with 
pion absorption and scattering on nuclei |18|, |19| , while the form [g§ + K 2 ]/[q 2 cm + k 2 } 



takes into account the A off-shell correction and has been used for calculating vac- 
A + N cross sections 



uum N + N 
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and is also frequently included in vacuum 
A widths incorporated in transport simulations [17], |2l|. Since the two form factor 
choices, FF1 and FF2, lead to somewhat different in-medium properties, we have 
chosen to test both forms. 

At the NpN and NpA vertices interactions corresponding to those used at the 
pion vertices are used, and at baryon-hole vertices effective short-range interactions 
are included, moderated by the correlation parameters g' NN , g'j^^, and g' AA . 

A set of RPA equations for spin-isospin modes, corresponding to an infinite itera- 
tion of non-interacting pion (p- meson), nucleon-hole, and A-hole states, were derived 
in Ref. [0. From these equations eigenvectors and eigenenergies are obtained for the 
different spin-isospin modes. We also calculate total and partial A widths, iV7r — > A 
cross sections, iViV <-> AiV cross sections and A spectral functions within the RPA 
approximation (explicit expressions are given in Ref. [|12]j ). 

The results of the AiV -1 model depend on the model parameters which cannot 



be determined uniquely by comparison to experimental data. In Ref. [12| was chosen 
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FF1 & FF2 


FF1 & FF2 


FF1 


FF2 


m N c 2 = 0.940 
m A c = 1.23 
m^c 2 = 0.14 
m pC 2 = 0.77 


fNN = 1-0 

/at a = 2-2 
/Xa = 
A*" = 1.0 


Inn = 6-2 
/&a = 10.5 
/aa = 
A p = 1.5 


9'nn = 0-9 
^a = 0-38 
g' AA = 0.50 
A, = 1.5 


9'nn = 0-9 
9na = 0-43 
9aa = 0-40 
A 9 = 1.5 


= m N 


po = 0.153 fm" a 
T = 10~ 4 


Va-Vn = 0.025 p/po , p < p 
V a -Vn = 0.025 , p > p 



Table 1: Parameter values used in microscopic calculations (energy units: GeV). 



a particular set that could reproduce vacuum pp — > A ++ n cross sections and the 



imaginary part of the empirical A-nucleus spreading potential of Ref. |[22|| . Unfor- 
tunately, the established transport treatment that we have chosen to compare our 
medium-modified simulations to [17j does not contain any effective nucleon mass, 



as was utilized in the microscopic treatment of Ref. ||12|| . For this reason we have 



here ignored the effective mass in the microscopic calculations (i.e. we have taken 
m* N = mjy). In addition we are using different sets of form factors. To still reproduce 
the vacuum pp — > A ++ n cross section and the empirical spreading potential of Ref. 



22] as well as possible we have readjusted the g' parameters (while keeping other 
parameters as in Ref. [0). With the parameter choice in Table |IJ, we find that FF1 
over predicts Im V^ prea d somewhat, while FF2 falls below the empirical points. Fur- 
thermore, FF1 yields a slightly larger pp — > A ++ n vacuum cross section than FF2. 
Therefore, we will often show results obtained with both form factor choices. 

The microscopic calculations are performed in a system where the medium is at 
rest and the obtained dispersion relations and decay widths refers to this medium 
frame. When incorporating the in-medium quantities into the transport formalism, 
this can be effectively taken into account of by defining a local medium frame as the 
frame where the local flow velocity vanishes. The pionic HamiltoniansQ are deduced 
from the density-dependent dispersion relations in infinite nuclear matter by a local- 
density approximation and are parameterized and utilized in the local medium frame. 

The pionic modes are created in the A decays which are governed by partial A 
widths. As illustrated in Fig. [1J, the in-medium widths are different than the free 
width. This is mainly because of different phase space and because the pion com- 
ponent on the pionic modes no longer is unity; there are also AA^ 1 and NN~ X 
components. Note particularly that the width for the decay to the lower mode de- 
creases and approaches zero for invariant masses above m ~ 1.3 GeV/c 2 . This is 
because the pion component vanishes and the collectivity disappears when the mode 
enters the band of non- collective AA^ _1 modes. This differs from the findings in Ref. 

TOU which did not include a AA^ 1 continuum. Otherwise Ref. |lfj] found qualitatively 



similar partial widths, with quantitative differences depending on the specific form 
factors and parameter sets. Another feature of the A width in the medium is that it 



1 In this work we explicitly propagate both the lower and the upper pionic modes, as was also done 
in Ref. JlTJ]. This differs from Refs. |llj where only a single average pionic mode was incorporated. 
The importance of treating the two pionic modes equally was discussed in Refs. ^] . 



3 




Figure 1: Partial A decay widths for the lower pionic mode tti (upper panel) for zero 
A momentum, and pion absorption cross sections for the process pirf — > A ++ (lower 
panel). The cross sections are obtained for the special case when the pnf cm. frame 
is equal to the rest frame of the medium. The curve labeled Kitazoe represents the 
parameterization in Ref. [21], and the empirical cross section for 7r + p scattering are 



estimated from Fig. 2.2 in Ref. 24 



depends explicitly on the A momentum in the medium, in addition to the A energy. 
This is especially important for the partial widths. To this end we calculate and 
store in a large table the partial A decay widths for discrete values of invariant mass, 
momentum, density and angle. The A then decays in the medium frame according 
to linear interpolation in the table of decay widths (for details see Ref. [pf]). 

The pionic modes may be reabsorbed through the process txN —>■ A. In the stan- 
dard transport treatment, the vacuum reabsorption cross section may be written as 
an interaction factor times a Breit-Wigner factor containing the free A width. The 
in-medium pion absorption cross section has a similar structure. For the transport 
treatment we make use of the fact that the interaction factor is also used for calculat- 
ing the partial A width, and we use the full in-medium A width in the Breit-Wigner 
factor. Examples of Nn — > A cross sections for FF1 are presented in Fig. [|. 

The transport simulations are based on the hadronic transport model of Li and 



Bauer [17]. In this model the dynamical evolution of a heavy-ion collision is described 



by a set of coupled transport equations, explicitly treating nucleons, A isobars and 
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pions. In this treatment the inelastic NN cross section is parameterized following 



Ref . |25[ . In the medium- modified simulations we will replace the VerWest and Arndt 
vacuum cross sections for the processes NN «-> AiV by the density dependent cross 
sections obtained from the microscopic calculations. These calculated cross sections 
constitute in vacuum a somewhat cruder approximation than the fit obtained in Ref. 



25] where three independent functions 0777 were used. Since our aim in this work is to 



test in-medium effects, we will compare to transport simulations based on the model 



described in |F7| , but with the VerWest and Arndt NN <-» AiV cross sections replaced 



by our calculated cross sections at zero density. This we will refer to as a standard 
simulation. The microscopic cross sections are calculated for the case when the cm. 
system of the colliding particles coincides with the rest system of the medium. To 
keep the transport treatment as simple as possible these cross sections are then used 
also when this cm. system differs from the local medium frame. 

One important medium modification of the NN — > AiV cross section originates 
from the A spectral function which is modified in the medium. This spectral function 
contains the full in-medium A width, which is calculated self consistently. The A 
spectral function is also used in the transport treatment to determine the invariant 
A mass, once a NN — ► AiV collision has occured. The total in-medium A width 
leads at high densities to a substantial reduction of the NN — > AiV cross section for 
moderate and large cm. energies and to an enhancement for low cm. energies. This 
differs from the results found in Refs. [[8], |2(| where an enhancement was found for 



large densities. However, as noted in Ref. ||, the result is very sensitive to the g' 
parameter values. In Ref. [T2(J it was shown that the enhancement found in Ref. || is 
weakend by the inclusion of the full A width in the pion polarization function. Since 
our calculations employ large values of the g' parameters^ and in several aspects differ 
from Refs. ||, J2fj, we find a net surpression of a(NN — > AiV). Most importanly, we 
include selfconsistently the full in-medium A width (including the A — > N + NN" 1 
contribution which is important at low A mass and high densities) and we perform 
an integration over A masses including the in-medium A spectral function. 

The cross section for the process iVA — > NN is obtained and implemented anal- 
ogously to the NN — > AN cross section. Note though that the absence of the A 
spectral function simplifies the treatment. However, another complication arises be- 
cause the cross section also depends on the invariant mass of the colliding A, in 
addition to the cm. energy. We have empirically found that a(y/s; m) « a{\fs i ] ttia), 
with y/s' = y/s — 0.778(m — jjia) yields a good approximation. The iVA — > NN cross 
sections show only a weak dependence on the nuclear density. 



3 Results 

In this section we present results from BUU simulations containing the in-medium 
quantities as discussed in Sect. §. The purpose of these simulations is to elucidate the 
effects of the included in-medium properties. Therefore, to make the treatment simple 
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Since the g' parameters are multiplied by form factors, there is not a one-to-one correspondance 



between the values used in Refs. |q, 26 and those used here 
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and the signals as clean as possible, we have ignored some processes that may be of 
equal importance when comparing to experiment. We have performed simulations 
for central 197 Au + 197 Au collisions at 1.0 A GeV and we have obtained results from 
standard simulations and from three different medium-modified simulations, denoted 
M1-M3. Ml contains in-medium modifications of only pion dispersion relations and 
partial A decay widths. M2 contains in addition to Ml also the in-medium pion 
reabsorption cross section a(Nn — > A). Finally, in M3 also the in-medium NN <->• 
AiV cross sections are added. All simulations were performed with 250 test particles 
per nucleon utilizing a mean field U(p) = -0.218(p/p ) + 0.164((p/p ) 4/3 GeV. 



1 A GeV Au+Au 




1 1.2 1.4 

Invariant A mass, m (GeV/c 2 ) 

Figure 2: Distribution of invariant mass for A's produced in NN collisions. 

The production rate of A isobars from NN collisions depends on the NN <-» AiV 
cross sections^ and the density and energy distributions of nucleons. When the density 
dependent in-medium NN — > AN cross sections are incorporated (simulation M3) 
the NN — ► AiV rates are reduced, due to a smaller in-medium cross section at high 
densities. For FF1 this effect is small, due to a strong enhancement of the cross 
section at low A/i-energies, while the reduction for FF2 is substantial. 

The in-medium NN — > AiV cross section, together with the A spectral function 
also lead to strong modifications of the invariant A mass distributions. Figure |2| 

3 The standard simulations with FF1 and FF2 yield slightly different results due to slightly 
different input vacuum NN <-> AiV cross sections. 



6 



shows that for simulation M3, there is a strong reduction of A isobars with invariant 
mass around = 1232 MeV, and an enhancement at lower masses. Note partic- 
ularly that the invariant mass distribution is non-vanishing also for masses below 
m = rriN + rn^. These A's cannot decay to a pionic mode and the only reabsorption 
channel is NA — > iVW. As a result, the net number of pions is somewhat reduced. 
The effects found with FF1 are similar to FF2 but show more pronounced and pecu- 
liar enhancement of low mass A's. This originates from the microscopic calculations 
which for FF1 show strong effects of pion condensation for densities around twice nor- 
mal nuclear density (while FF2 does not lead to pion condensation at the densities 
probed in the applied transport simulations). In the microscopic calculations, this 
phenomenon manifest itself as large enhancements of the A decay width to low-energy 
nucleon-hole modes which in turn leads to a strong enhancement of the A spectral 
function around m « 1 GeV and corresponding effects in the NN — > AiV cross sec- 
tion. These effects survive the transport treatment and we find strong enhancements 
of the number of produced A's in the invariant mass region around 1 GeV. The net 
effect on the number of produced pions is a slight reduction. 

Inspecting the rate of produced pions from A decay (upper panel of Fig. |]), we 
find that for FF1 there is in simulation Ml a slight reduction as compared to the 
standard treatment, due to a reduction in the partial A width at high densities^. 
FF2, which for m > ira has only about 60% of the corresponding width obtained 
for FF1, gives a strong reduction in the pion production rate. The further reduction 
in the A decay rate for simulations M2 to M3 reflects that there are fewer produced 
A's in the system. The simulation M2 illustrates this effect due to a reduction of the 
process Ntt — > A while for simulation M3 the reduction is due to fewer produced A's, 
in the mass region m > + m n , from NN collisions. 

The pion reabsorption rate (lower panel of Fig. |3|) depends on the density of pions 
and nucleons and on the reabsorption cross section. The reduction in the rate for 
simulation Ml is due to the fact that fewer pions are produced from A decay. There 
is a further strong reduction when including also the in-medium Nn — ► A cross 
section (M2). This is because this cross section is suppressed at high densities (see 
Fig. |l|). The pion reabsorption rate is reduced even further for simulation M3. This 
is a secondary effect of the fewer produced A's which then yield fewer pions. 

For FF2 the net number of A isobars and pions present in the system at different 
times varies much less with the in-medium input than the production and absorption 
rates. This is because different effects balance each other and because the production 
and absorption depends on particle densities which yields an "automatic" regulation 
of the net number. Note particularly that the enhancement of pions due to the low 
reabsorption cross section at high densities is compensated by a low production of A 
isobars by a likewise low NN — > AiV cross section at high densities. This compen- 
sation is not accidental as it depends mainly on the A in-medium spectral function 
which is included in both cross sections. The net effect is a reduction of the number of 
emitted pions of only about 3-10% for simulation M3 as compared to the result from 

4 For the medium-modified simulations we find that about 95% of the emitted pions originate 
from the lower pionic mode. 
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Figure 3: The time dependence of the rates for A decay (upper panel), and the time 
dependence of the rates for pion reabsorption (lower panel). 



the standard simulation. For FF1 the number of present A's varies somewhat more 
with the different types of medium-modified simulations. These variations originates 
from the microscopic onset of pion condensation, as briefly discussed above. 

In-medium effects show up also in the pion energy spectrum, although to a rather 
small degree. Transverse momentum spectra da/prdpr have been obtained at impact 
parameter b = 1.0 fm assuming that all collisions with an impact parameter up to 
&max=2r A^ contribute equally. The effect seen in the simulations incorporating 
the in-medium properties, as compared to the standard simulations (see Fig. [|), is a 
modest but significant enhancement at low transverse momenta (pt < 250 MeV/c) 
and a reduction at higher momenta (pr >350 MeV/c), corresponding to a reduction 
of the effective transverse temperature. We find the enhancement of low-energy pions 
encouraging, even though this effect cannot alone explain the enhancement seen in 
experimental spectra (Fig. 1 of Ref. |27|]). It should be emphasized that the obtained 
spectra are not entirely suitable for quantitative comparison with experimental spec- 
tra partly because of the simple impact-parameter averaging employed, and partly 
because the incorporation of higher nucleon resonances and additional medium effects 
might be of importance. Note, for example, that higher-lying resonances contribute to 
both low and high energy pions through the two-pion and one-pion emission channels, 
respectively f28fl. 
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Figure 4: The transverse momentum spectrum for neutral pions in the central 
rapidity interval —0.20 < y cm < 0.20. The error bars represent the statistical errors. 
The results obtained for FF2 are very similar. 

4 Summary 

In this letter, we have presented selected results from transport simulations incorpo- 
rating in-medium properties of pions and A isobars. In particular, we have included 
in-medium pion dispersion relations, partial A widths, pion reabsorption cross sec- 
tions, NN <-> AiV cross sections and A spectral functions. These in-medium quanti- 
ties have been calculated microscopically from the AiV -1 model presented in Ref. |12| 



and incorporated into the transport formalism of Li and Bauer [[17]] by a local-density 
approximation. The medium-modified simulations presented in this work show strong 
effects on properties not directly observable during the collision process, such as pion 
and A production and reabsorption rates, but only minor effects on spectra of emitted 
pions. This is rather reasonable since most of the emitted pions are produced at the 
surface at low densities where the in-medium effects are quite small. However, in an 
energetic nucleus-nucleus collision also other particles, not studied in this work, are 
produced in multistep processes, where A's and pions act as intermediate particles. 
Thus, modified properties of A's and pions in the nuclear medium might be important 
to consider when studying production of such secondary particles. 
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